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ABSTRACT 

We derive fundamental parameters for the old, rich open cluster NGC 1245 using BVI photometry 
from the MDM 1.3m and 2.4m telescopes. Based on detailed isochrone fitting, we find NGC 1245 
has a slightly sub-solar metallicity, [Fe/H] = —0.05 ± 0.03 (statistical) ± 0.08 (systematic), and an 
age of 1.04 ± 0.02 ± 0.09 Gyr. In contrast to previous studies, we find no evidence for significant 
differential reddening. We determine an extinction of Ay = 0.68 ± 0.02 ± 0.09 mag and a distance 
modulus of (m — M)o — 12.27±0.02±0.12 mag, which corresponds to a distance of 2.8±0.2 Kpc. We 
derive a logarithmic mass-function slope for the cluster of a = —3.12 ± 0.27, where a Salpeter slope 
is a = —1.35. Fits to the radial surface-density profile yield a core radius of r c = 3.10 ± 0.52 arcmin 
(2.57 ± 0.47 pc). NGC 1245 is highly relaxed and contains a strongly mass segregated population. 
The mass function for the inner cluster has a very shallow slope, b — —0.56 ± 0.28. In contrast, the 
outer periphery of the cluster is enriched with low mass members and devoid of high mass members 
out to the tidal radius, r t = 20 arcmin (16.5 pc). Based on the observed surface-density profile and 
an extrapolated mass function, we derive a total cluster mass, M = 1300 ± 90 ± 170M Q . 

Subject headings: stars: fundamental parameters — stars: luminosity function, mass function — open 
clusters and associations: individual (NGC 1245) 



1. INTRODUCTION 

Open clusters are excellent laboratories for many dif- 
ferent aspects of astrophysics. First, open clusters form 
a coeval set of stars with homogeneous properties. This 
homogeneity makes it possible to determine with rela- 
tive e ase the age and metallicity of the cluster (lYi et al. I 
HHI). Second, owing to their relatively small relax- 
ation times, open clusters provide an opportunity to 
study s tellar systems in various st ages of dynamical evo- 
lution ijBinnev fc Tremaind 11987)) . Finally, theory pre- 
dicts the star formation process within an open cluster is 
strongly affected by dynamical interactions, supernova e 
explosions, and UV radiation l)Adams fe M vers"! 12001). 
Characterizing the fraction of cluster members with com- 
panions from stellar to planetary masses for a cluster in 
comparison to the field provides valuable insight into how 
these processes affect stellar and planetary formation. 

The large angular size (tens of arcminutes to several 
degrees) of open clusters has traditionally made it diffi- 
cult and time consuming to observe a substantial fraction 
of an open cluster's members. However, recently avail- 
able large-format CCD imagers now allow for complete 
studies of an open cluster with relatively small amounts 
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of observing time. In this paper we present wide-field 
BVI photometry of the open cluster NGC 1245 using the 
MDM 8K Mosaic imager on the MDM 2.4m Hiltner tele- 
scope and supplemental BVI photometry obtained un- 
der photometric conditions on the MDM 1.3m McGraw- 
Hill telescope. 
NGC 1245 is a rich, old open cluster with 

an approximately so lar me tallicity po p ulatio n 

Panes. Tillev. fc Lvmral Il988l: IWee fc Lee I EHf). 
The cluster is located toward the galactic anticenter 
£ = 147°, b = -9°) at a distance, R ~ 2 .5 kpc 
Janes. Tillev. fc Lvnera lll98^: IWee fc Leelll996fl . With 
its relatively large Galactocentric distance, NGC 1245 
is particularly usefu l in constr a ining the Galacti c 
metallicity gradient <)Friel I Il995t IWee fc Lee I Il996|) . 
Here, we improve on the observations for this cluster 
by covering a six-times greater area, obtaining the first 
CCD /-band photometric data, and acquiring the first 
CCD data under photometric conditions. With these 
significant improvements over earlier studies, we are able 
to determine more precisely the physical parameters of 
the cluster as a whole, rather than just the inner cluster 
core. In addition to the observational improvements, 
we improve the interpretation of the observations by 
calculating the physical parameters of NGC 1245 using 
a chi-square fit to the observations, and we make every 
effort to quantify the systematic errors in the derived 
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physical parameters. 

Because NGC 1245 is a rich open cluster and read- 
ily observable from our primary observing site, it is the 
first target of the Search for Transiting Ex trasolar Plan- 
ets in Stellar Systems (STEPSS) project l|Burke et al. I 
12003ft . The goal of the STEPSS project is to assess the 
frequency of close-in extrasolar planets around main se- 
quence stars in several open clusters. Part of the BVI 
data we are presenting in this paper were taken during a 
19 night observing run to search for transiting extrasolar 
planets in NGC 1245. This transit data consist of /-band 
only observations obtained every seven minutes with the 
MDM 8K Mosaic imager on the MDM 2.4m telescope. 
Results of the search for extrasolar planets in NGC 1245 
using this transit data will be presented in an upcoming 
paper (B. S. Gaudi et al., in preparation). An additional 
paper presenting the variable star content of our NGC 
1245 field based on light curves from the transit data is 
also in progress (J. Pepper et al., in preparation). The 
main goal of this paper is to determine the fundamen- 
tal physical parameters of NGC 1245. With the known 
properties for this cluster and subsequent clusters in the 
STEPSS project, we will gain insight into how metallic- 
ity, age, and stellar density affect planetary formation, 
migration, and survival. 

We describe our observations and data reduction in 
33 In 33 w e derive the physical parameters of NGC 
1245 from isochrone fits to the BVI photometry. We 
derive the radial profile of NGC 1245 in <21 the mass 
function in 33 and the total mass in 33 We summarize 
and conclude in 33 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

We observed NGC 1245 on two occasions. The first set 
of observations was obtained in November 2001 using the 
MDM 8K mosaic imager on the MDM 2.4m Hiltner tele- 
scope. The MDM 8K imager consists of a 4x2 array o f 
thinned 2048x4096 SITe ST002A CCDs l)Crotts I [200 11. 
This instrumental setup yields a 26'x26' field of view and 
a 0.36" per pixel resolution employing the 2x2 pixel bin- 
ning mode. Table 2] shows for each night of observa- 
tions the number of exposures in each filter, exposure 
time, median full width at half maximum (FWHM), and 
a brief comment on the observing conditions. This ob- 
serving run is optimized for a transiting extrasolar planet 
search, thus there are 960 I-band images obtained during 
this run that are not listed in Tabled Since the empha- 
sis of the present study is the color-magnitude diagram 
of NGC 1245, we chose only the I-band images on each 
night that were taken concurrent with the other two fil- 
ters and the same number of images to have similar signal 
to noise characteristics for each filter. None of the nights 
were photometric. Therefore, we reobserved NGC 1245 
in February 2002 using the MDM 1.3m McGraw-Hill tele- 
scope with the 2048x2048 "Echelle" imaging CCD. This 
CCD provides a 17'xl7' field of view with 0.5" per pixel 
resolution . Several night s were photometric, and we use 
Landolt ijLandolt 111992ft standard star observations to 
calibrate our photometry. Table |3 details the images 
taken with the MDM 1.3m telescope. The relative place- 
ment between the 1.3m and 2.4m field of views is shown 
in Figure n 

Figure ^ shows a false-color image of NGC 1245 and 



shows the overlap between the MDM 2.4m and 1.3m 
datasets. The grid of eight numbered rectangles shows 
the field of view for the eight CCDs that make up the 
MDM 8K Mosaic imager. The large inscribed square is 
the field of view of the MDM 1.3m observations. The 
blue, green, and red color channels for the false-color im- 
age result from a median combined image of the BVI 
passbands, respectively. The circle denotes the cluster 
center derived in Section 0] 

2.2. Data Reduction 

We use the IRAF 1 CCDPROC task for all CCD pro- 
cessing. We describe the 2.4m data reduction first. The 
stability of the zero-second image over the course of the 
19 nights allows us to median combine 95 images to deter- 
mine a master zero-second calibration image. For mas- 
ter flat fields, we median combine 11, 19, and 66 twilight 
sky flats taken throughout the observing run in the BVI 
passbands, respectively. We quantify the errors in the 
master flat field by examining the night to night variabil- 
ity between individual flat fields. The small-scale, pixel- 
to-pixel variations in the master flat fields are ~ 1%, and 
the large-scale, illumination-pattern variations reach the 
3% level. The large illumination-pattern error results 
from a strong sensitivity in the illumination pattern to 
telescope focus. However, such large-scale variations do 
not affect photometry where the point-spread-function 
(PSF) scale is the most relevant. 

We use the DAOPHOT package ijStetson 11987ft within 
IRAF for PSF fitting photometry, and we perform the 
photometry on each frame individually. To identify 
stars and to calculate stellar positions, we designate 
the highest signal-to-noise /-band image as input to the 
DAOFIND task. We transform the stellar positions de- 
termined on the high signal-to-noise /-band image to all 
the other images as initial guesses for photometry. To 
calculate the model PSF, we select 25 bright, isolated 
stars evenly distributed on the reference image. We use 
the same PSF stars for all images. Slight offsets be- 
tween images result in at least twenty stars in the PSF 
model calculation for all the images. Modeling the spa- 
tially constant PSF model involves three iterations with 
faint neighbors to the PSF stars being subtracted from 
the image between iterations. Using the resulting PSF 
model, the IRAF ALLSTAR task calculates the instru- 
mental magnitudes by fitting the PSF model to all stars 
simultaneously. The reference image photometry defines 
our instrumental magnitude system, and we apply a zero- 
point offset to the individual image photometry to align 
them with the reference image photometry. A robust 
weighted mean using the errors output by ALLSTAR 
gives the final instrumental magnitude. 

In order to place the MDM 2.4m instrumental mag- 
nitudes onto a standard system, we use the IRAF FIT- 
PARAMS task to determine the transformation to the 
photometric MDM 1.3m dataset. The rms scatter for an 
individual star around the best fit photometric transfor- 
mation is ^0.03 mag in all passbands. This error repre- 
sents how well an individual star approximates the pho- 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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tometric system defined by the 1 .3m data described be- 
low. The data taken on CCDs 1 and 3 of the MDM 8K 
mosaic imager suffered from a nonlinearity that reaches 
a level of 0.14 mag near full well. With the aid of the 
1.3m dataset, this nonlinearity is correctable such that 
the errors in the photometric calibration for these CCDs 
are similar to the other CCDs not affected by the nonlin- 
earity. The bottom panels of Figure show the resulting 
color magnitude diagrams (CMDs) for the MDM 2.4m 
dataset. The stars shown have sample standard devia- 
tions in the BVI photometry of less than 0.08 mag and 
at least 3 (B), 3 (V), and 5 (I) photometric measurements 
return from the photometry package without numerical 
convergence errors. The error bars in Figure |21 show the 
median of the sample standard deviations in the photo- 
metric color determination for stars in bins of 0.5 mag 
at V. These error bars represent the expected deviation 
from the mean of a color measurement from a single set 
of BVI exposures with similar signal to noise properties 
as our dataset. Our ability to determine the mean color 
for a star requires dividing these errors by the yf~N of the 
number of measurements. 

An alternative method to analyze the MDM 2.4m 
dataset is to combine the images and perform the pho- 
tometry on this master combined image. The combined 
image enhances the signal to noise of the faintest stars. 
This enhanced signal to noise allows DAOFIND to de- 
tect and ALLSTAR to fit the model PSF for significantly 
fainter stars at the expense of having to trust the error 
output from the photometry package. We median com- 
bine 15, 30, and 31 BVI images, respectively, and per- 
form the identical photometry procedure outlined above. 
Initial guesses for the stellar positions come from the 
combined I-band image. The instrumental magnitudes 
are fit to the calibrated MDM 1.3m dataset with a sim- 
ple zeropoint and color term. The resulting CMDs are 
shown in Figure As one can see these CMDs go signif- 
icantly fainter, but the field star contamination becomes 
increasingly a problem, especially for the V,B — V CMD. 
The low contrast of the cluster's main sequence against 
the field star contamination limits the utility of these 
deeper CMDs in determining the physical parameters of 
NGC 1245. In this study we only use these deep CMDs 
for extending our mass function toward fainter magni- 
tudes (see Section |SJ| . 

The analysis procedure for the MDM 1.3m dataset is 
similar to the MDM 2.4m dataset. The master zero- 
second image consists of a median combination of 49 im- 
ages, and a master sky flat consists of a median com- 
bination of 12, 15, and 4 images for BVI passbands, 
respectively. The PSF has strong spatial variations, so 
we choose the same set of ~ 80 stars evenly distributed 
across a reference image to calculate a quadratically vary- 
ing PSF model for each image individuall y. We employ 
the s tand-alone ALLFRAME program ijStetson et alJ 
1998) to calculate the instrumental magnitudes for all 
images simultaneously. The aperture corrections also are 
spatially variable. We fit a third order 2-D polynomial 
to the aperture corrections calculated on the evenly dis- 
tributed PSF stars. The IRAF XYZTOIM task provides 
the best fit solution for the 2-D polynomial aperture cor- 
rection model. The standard deviation of the aperture 
correction is ~ 0.02 mag. 



Aperture photometry m easurements of the Landolt 
field SA101 ^Landolt l ll9921 taken at three different air- 
masses, X =1.2, 1.5, and 1.9, on two separate photomet- 
ric nights (see Table 0) provide a determination of the 
airmass coefficient. The airmass coefficients for the two 
nights are statistically identical; thus, we average both 
airmass coefficient determinations for the final calibra- 
tion. Fixing the airmass coefficients, we calculate color 
terms for each night using observations of the Landol t 
fields G44-27, G163-50/51, and PG1407 {Landolt 11992ft . 
The color coefficients from the individual nights are also 
statistically identical, and we use the mean of both color 
coefficient determinations for the final calibration. The 
standard deviations in the calibrations are 0.015, 0.015, 
and 0.02 mag for the BVI passbands, respectively. This 
represents how well an individu al star approxim ates the 
photometric system outlined bv lLandolt 1 1)1992|) . 

There are 2, 2, and 3 images in the BVI passbands, 
respectively, of NGC 1245 that are photometric. After 
applying the aperture corrections to the ALLFRAME in- 
strumental magnitudes for the photometric data, we take 
a weighted mean of the measurements and then apply the 
photometric calibration to derive an initial standard-star 
catalog. The final standard-star catalog is obtained by 
applying the color term and aperture corrections to all 
nonphotometric MDM 1.3m data, fitting for the zero- 
point offset to the initial standard-star catalog, and cal- 
culating the robust weighted mean of all measurements. 

In the process of determining the zero-point fit for the 
nonphotometric MDM 1.3m data to the initial standard- 
star catalog, we uncovered a 1.5% amplitude residual 
that correlates with spatial position on the CCD. The 
residual is negligible over most of the CCD, but the resid- 
ual abruptly appears for stars with a position x < 500 
in pixel coordinates. The residual linearly grows toward 
decreasing x position on the detector reaching the 1.5% 
level at the edge. Due to the abrupt onset of this resid- 
ual, even a cubicly varying PSF model did not reduce 
the amplitude of the residual. To reduce the effect of this 
systematic spatial variability, we limit the zero-point fit 
for the nonphotometric MDM 1.3m data to stars in the 
middle of the detector defined by 500 < x < 1500 and 
500 < y < 1500 in pixel coordinates. 

The top panels in Figure El show the resulting CMDs 
for the MDM 1.3m dataset. Only stars that have at least 
3 photometric measurements return from the photome- 
try package without numerical convergence errors and a 
sample standard deviation < 0.08 mag in all three pass- 
bands are shown. 

2.3. Comparison with Published Photometry 

There is some confusion in the litera ture regarding the 
photo metric calibration of NGC 1245. iCarraro & Patatl 
(1994) perform the first CCD detector study of NGC 
1245. They calibrate their photometry by compari- 
son w ith the photoelectric observatio ns of [Hga^^^ 
(1961). Unfortunately, the data of Carraro & Patat 
(1994) available on the widely used WEBDA open clus- 
ter database 2 is incorrect since it lacks the photomet- 
ric zeropoint as used in publication (G. Carraro pri- 
vate communic ation). Thu s, usin g the incorrect data 
from WEBDA, iWee fc Lee I l(T99^ concluded the cali- 

2 http://obswww.unige.ch/webda/ 
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bration of ICarraro fc Patat I lfl994|) was incorrect. Our 
independent photometric calibration agrees with the ear- 
lier photoelectric studies. For 5 and 3 stars i n com mon 
with the photoelec t ric stu dies of lHoag et al. I (^1961^1 and 
I.Iennens fc Helferl l)1975lk respectively, we find V-band 
differences of 0.03 ± 0.06 and -0.012 ± 0.007 mag, re- 
spectively, and (B — V) color differences of —0.01 ± 0.01 
and —0.001 ±0.02 mag, respectively. The (B — V) colors 
of the stars used in the comparison span the range from 
0.6-1.2 mag. Thus, the photometric calibration for NGC 
1245 is well determined . 

The recent study by ISubrama niaml 1 20031) does not 
find any differences with the ICarraro fc Patatl l)1994f ) 
photometry in the V band but does find a linear trend in 
the (B — V) color difference as a function of V magnitude 
with an amplitu de of 0.35 mag. It is unc lear if the com- 
parison with th e lCarraro fc Patatl l)1994l) photometry in 
ISubramaniam I l)2003j) is with the incorrect data from 
WEBDA. However, we do not find any (B-V) c olor differ- 
ence trends with the lCarraro fc Pa tat I lll99 4|) dat a. Ad- 
ditionally, inspection of the ISubramaniam I l|2003j) CMD 
finds the giant clump is 0.35 mag redder in (B — V) than 
in our CMD. The above comparisons call into question 
the ISubramaniaml l)2003j) color calibration. 

3. PHYSICAL PARAMETERS 

NGC 1245 is located in the Galactic plane (b = -9°), 
and therefore, the observed field contains a large num- 
ber of foreground and background disk stars compli- 
cating the process of fitting theoretical isochrones to 
the CMDs. We attempted to statistically subtract 
the field-star population with a pro cedure similar to 
iMiehell. Saraiedini. fc French I l)1998j) . Unfortunately, 
because the contrast of the main sequence above the 
background fiel d stars was relatively low (in compariso n 
to the work by [Mighell, Saraie dini. fc French (1998)), 
we were unable to achieve a satisfactory background field 
subtraction. 

Rather than statistically subtracting field contamina- 
tion, we select stars for the main sequence isochrone fit 
by the following procedure. Since binaries and differen- 
tial reddening affect the (B — V) color of a star less than 
its (V - 1) color, the V, (B - V) CMD has a tighter main 
sequence than the V, (V — I) CMD. The tighter main 
sequence in the V, (B — V) CMD has a higher contrast 
above the contaminating field stars. Thus, to select main 
sequence cluster members for the isochrone fit, we trace 
by eye the blue main sequence edge in the V, (B — V) 
CMD. Shifting this trace redward by 0.12 mag defines 
the red edge selection boundary. Finally, shifting the 
trace 0.02 mag blueward defines the blue edge selection 
boundary. These shifts were selected by eye to contain a 
high fraction of the single-star main sequence. The ver- 
tical light solid lines in Figure 0] show the resulting main 
sequence selection boundaries in the V, (B — V) CMD. 

Since the difference in color between the red-giant 
clump and the main sequence turnoff in the CMD places 
strong constraints on the cluster age, we designate a 
box in the V, (B — V) CMD to select red giant mem- 
bers for the isochrone fit. The region with boundaries 
13.7 < V < 14.4 and 1.09 < (B - V) < 1.21 defines the 
red-giant clump selection. The light solid box around the 
red-giant clump in the V, (B — V) CMD in Figure|3|shows 
the selection criteria for these stars. The red-giant clump 



is saturated in the MDM 2.4m dataset, and therefore, we 
fit isochrones to the MDM 1.3m photometry only. There 
are 1004 stars that meet our main sequence and red-giant 
clump selection criteria. 

For the isochr one fitting, we use the Yale-Yonsei 
(Y 2 ) isochrones (lYi et al. I 120011) . which employ the 
iLeieune. Cuisinier. fc Buser I <|1998i) color calibration. 
The Y 2 isochrones provide an interpolation scheme to 
calculate isochrones for an arbitrary age and metallic- 
ity within their grid of calculations. For a given set of 
isochrone parameters (metallicity, age, distance modu- 
lus, Ay, and -Ry), we define the goodness of fit as 

Xtot = ^2xi, (1) 

i 

where the sum is over all stars selected for isochrone fit- 
ting, and xf is the contribution from star i, 

Xi = (Bpred(mi) — B b s ,i) 2 + (V pre d(mi) — V bs,i) 2 + (Ipred(mi) — Iobs,i) 2 

In the preceding equation, BVI pre d are the magnitudes 
predicted by the isochrone, using the stellar mass, m„ as 
the independent variable, and BVl b s ,i are the observed 
stella r magnitudes. We use the Brent minimization rou- 
tine l)Press et al. I l200lt) to determine the stellar mass 
that minimizes xf for each star. 

Theoretical isochrones and color calibrations con- 
tain systematic uncertainties that make it impossible 
to find a consistent fit to CMDs of the best stud- 
ied open cl usters over a wide range of colors and 
magnitudes (Ide Bruiine. Hoogerwerf. fc de Zeeuw 1120011: 
iGrocholski fc SanHechni I l2003|h In order to avoid 
overemphasizing the main sequence turnoff, where the 
photometric errors are the smallest, we adopt an equal 
weighting to all stars. 

To find the best fit isochrone parameters, we per- 
form a grid search over age and metallicity. At fixed 
age and metallicity, we us e the Powell mult idimcn- 
sional minimization routine l|Press et al. I l200"ll) to de- 
termine the best fit distance modulus and reddening, 
Ay, by minimizing Xtot- To determine the reddening 
in the other passbands, As and Ai, we must assume a 
value for Ry, the ratio of total to selective extinction 
(jCardelli. Clavton. fc Mathis 111989)) . For our best fit so- 
lution, we fix Ry — 3.2. We discuss the reasoning for this 
choice and its impact on the results in the Section l3~Tl 
We calculate xiot over a 40x40 evenly-log-spaced grid of 
metallicity and age that ranges —0.26 < [Fe/H] < 0.13, 
and 0.7943 < Age(Gyr) < 1.4125. 

We determine our confidence limits on the model pa- 
rameters by scaling the resulting chi-square statistic, 

2 2 

A 2 Xtot Xrain /n\ 

X ~ v 2 . Iv ' 1 ' 

where X m % n is the minimum xiot an d v is the num- 
ber of degrees of freedom. Figure [S] is a contour plot 
showing the confidence region for the joint variation in 
metallicity and age. The three solid contours represent 
the 1-, 2-, and 3-cr confidence limits corresponding to 
Ax 2 = 2.3,6.14, and 14.0 for the two degrees of free- 
dom, respectively. To further refine the best fit isochrone 
parameters, we ran another 40x40 metallicity and age 
grid at finer resolution: —0.17 < [Fc/H] < 0.03 and 
0.9550 < Age(Gyr) < 1.2303. We quantify the 1-a error 
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bounds by fitting a paraboloid to the A\ 2 surface as a 
function of the four fit parameters [Fe/H], Age, (m— M)o, 
and Ay. The 1-er errors are then the projection of the 
Ax 2 = 1.0 extent of this paraboloid on the parameter 
axes. We note that this procedure assumes that the A% 2 
surface near the minimum is parabolic, which we find to 
be a good approximation. Also, because we minimize on 
Ay and (m — M)q at fixed age and metallicity, we are 
ignoring the additional variance from these parameters. 
Therefore, our errors are slightly underestimated. The 
resulting best fit parameters are [Fe/H] = —0.05 ± 0.03, 
Age = 1.035 ±0.022 Gyr, (m-M) = 12.27±0.02 mag, 
and Ay = 0.68 ± 0.02 mag. The best fit isochrone is 
overplotted on the V, (B - V) and V, {V - I) CMDs in 
Figure^ and the best fit metallicity and age are shown 
as the filled square in Figure 

The isochrone fits reasonably well in the V, (B — V) 
CMD, deviating slightly too red toward the lower main 
sequence. We investigate whether an adjustment to the 
photometric calibration color term can reconcile the ob- 
served main sequence and the isochrone by first calculat- 
ing the color residual between the stars and the isochrone 
at fixed V mag. A line fit to the color residuals as a func- 
tion of (B - V) color has a slope, b = -0.065 ± 0.007. 
This calculation assumes that the V-band color term is 
correct, and only the B-band color term needs adjusting. 
The error in determining the color term, 1 — a ~ 0.004 
in all passbands, suggests the large, b = —0.065, color 
term adjustment required to reconcile the observed main 
sequence and theoretical isochrone does not result from 
an incorrect determination of the color term alone. 

The fit is poorer in the V, (V - I) CMD. Visually, the 
best fit isochrone exhibits a tendency to lie redward of 
the main sequence in the V, (V — I) CMD. This tendency 
toward the red partly results from the broader and more 
asymmetric main sequence in the V, (V — 7) CMD. Addi- 
tionally, the overly red isochrone results from systematic 
errors in the shape of the isochrone. The isochrone fits 
best near (V — I) ~ 0.9 and deviates redward of the main 
sequence toward fainter and brighter magnitudes. A line 
fit to the (V — I) color residual between the stars and 
isochrone at fixed V mag as a function of (V — I) color 
does not reveal any significant trend. Again, as with the 
V, (B — V) CMD, the theoretical isochrones are unable to 
fit the detailed shape of the cluster main sequence from 
the turnoff down to the lowest observed magnitudes. The 
isochrone deviates from the observed main sequence by 
as much as 0.06 mag in (V — I) color at fixed V mag. 

3.1. Systematic Errors 

There are four main contributions to the systematic 
uncertainties in the derived cluster parameters: uncer- 
tainty in Ry, photometric calibration errors, theoreti- 
cal isochrone errors, and binary star contamination. We 
address the effect of each of these uncertainties on the 
derived physical parameters in turn. 

We can, in principle, determine Ry by additionally 
minimizing Xtot with respect to Ry. However, in doing 
so we find a best fit value of Ry = 2.3 and a metal- 
licity of [Fe/H] = —0.26. Val ues of Ry this low are 
rare in our galaxy (however, see lCould. Stutz. fc FrogeTl 
(2001)), and metallicity measurements based on medium 
resolution spectroscopic indices for NGC 1245 rule out a 
metallicity this low (J. Marshall et al., in preparation). 



We also can use infrared photometry fro m the 2MASS 
All-Sky Data Release Point-Source Catalog ijCutri et al. I 
2003) 3 to place more robust constraints on Ry. As 
a side note, 2MASS photometry does not reach nearly 
as faint as our BVI photometry. Thus, we do not in- 
clude it directly into our best fit solution for the phys- 
ical parameters of NGC 1245. To aid in determining 
Ry, we constrain the sample of stars employed in the 
isochrone fit to stars with V < 16, which corresponds 
to K s ~ 14.5 mag for turnoff stars. At this limiting 
magnitude, the typical photometric error in the 2MASS 
photometry is a Ks < 0.07. We transform the 2MASS 
K s pas sband magnitudes to the if-band system de- 
fin ed bv jBessell fc Brett! lll988l) f t he sy stem employed in 
the lLeieune. Cuisinier. fc Buser I l(1998j) color calibration 
used by the Y 2 isoc hrones) by adding 0.04 mag to the 
2MASS photometry ijCutri et al. 112003)) . The photomet- 
ric transformation of the 2MASS photometry from K s to 
K does not require a color-term l|Cutri et al. 1120031) . Us- 
ing the BVIK magnitudes of the V < 16 turnoff stars, 
we derive a best fit to the isochrones with Ry ~ 3.2. 
However, the resulting best fit age is ~ 1.2 Gyr, some- 
what older than the age of ~ 1.0 Gyr derived from the 
deeper BVI data alone. Conversely, fixing the age at 
~ 1.0 Gyr in our isochrone fit to the turnoff stars, we find 
Ry ~ 3.4 is the best fit solution. Figure^shows the best 
fit isochrone to the BVI data systematically deviates 
from the observed main sequence at the turnoff. Since 
the 2MASS photometry limits our BVIK isochrone fit to 
within 1.5 mag of the turnoff, we suspect the difficulty in 
determining Ry with an age consistent with our best fit 
solution arises from systematic errors in the isochrones. 
We therefore adopt a fiducial value of Ry =3.2 for our 
best isochrone fit but explore the impact a systematic 
uncertainty in Ry of 0.2 has on the cluster parameters 
in the following paragraph. 

We quantify the effect of the systematic uncertainty 
in Ry on the cluster parameters by recalculating the 
isochrone fit assuming Ry = 3.0, the l-er lower limit 
in our derived value of Ry. The resulting best fit is 
[Fe/H] = -0.1, Age=1.036 Gyr, (m - M) = 12.29 mag, 
and Ay — 0.68 mag. The open triangle in Figurc[5]shows 
the resulting best fit metallicity and age for Ry = 3.0. 
Since the (B — V) color is more sensitive to metallicity 
variations than the (V—I) color, the color difference be- 
tween the main sequence in the V, (B — V) CMD and 
V, (V — I) CMD is a good indicator of metallicity. Un- 
fortunately, for a fixed Ay, variations in Ry also alter 
the the relative (B — V) and (V — I) color of the main 
sequence. This degeneracy between Ry and metallicity 
in determining the main sequence color limits our ability 
to determine the metallicity. 

To assess the systematic error associated with the pho- 
tometric calibration uncertainty, we refit the isochrones 
assuming the /-band photometry is fainter by 0.02 mag. 
The resulting best fit is [Fe/H] = 0.0, Age= 1.021 Gyr, 
(m - M) = 12.31, and Ay = 0.62. The /-band offset 
solution is shown as the starred point in Figure 

We attempt to investigate the systematic uncertainties 
in the isochrones by fitting the Padova group isochrones 
ijOirardi et al. 1120001) to the CMDs. We find that these 
isochrones do not match the shape of the observed main 

3 http: / /www. ipac.caltech.edu/2mass/releases/allsky/doc/cxplsup. html 
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sequence. Systematic differences between the isochrones 
and data can reach up to 0.08 mag in (B — V) and (V — I) 
colors. The mismatch in shape to the observed main se- 
quence prevents precise isochrone fits with the Padova 
isochrones, and we therefore cannot reliably quantify the 
systematic errors in the isochrones from these fits. How- 
ever, we find that the solar-metallicity, 1 Gyr isochrone of 
the Padova group is broadly consistent with the observed 
data. 

As an alternative method to quantify the systematic 
errors in the isochrones, we calculate the physical param- 
eters for the well studied Hyades open cluster using an 
identical procedure to the NGC 1245 data. The Hyades 
CMD data are selected for membership based on proper 
motion and eliminated of binary sta r contaminants based 
on sp ectroscopic observations (see IPinsonneault et al. I 
120031 for a discussion of the data and membership se- 
lection). Using a metallicity and age grid with resolu- 
tion A[Fe/H] = 0.01 dex and Alog(Age) = 0.006 dex, 
we find best fit isochrone parameters of [Fe/H] = 0.10, 
Age = 670 Myr, (m - M) n = 3.23 mag, and Ay = 0.03 
mag. The recent study bv lPaulson. Sneden. fc Cochran I 
(2003) and many previous studies determine a value of 
[Fe/H] = 0.13 ± 0.01 for the metallici ty of the Hyades 
using high resolution spectroscopy. iPerrvman et al.~l 
(1998) find an age of 625 ± 50 Myr using an independent 
set of isochrones. The Hipparcos distance to the Hyade s 
is (m - M) = 3.33 ± 0.01 mag l|Perrvman et al. Ill998j) . 
and the extinction for the Hyades is commonly cited as 
negligible. The difference in the physical parameters of 
the Hyades from our isochrone fit to the more accurate 
determinations of these parameters provides the relative 
systematic error in our isochrone fitting technique due to 
systematic errors in the isochrones. 

The final systematic error source we address is the error 
resulting from contamination due to unresolved binaries. 
Unresolved binaries tend to lie redward of the single-star 
main sequence, and this binary "reddening" has a larger 
effect on the (V — I) color than the (B — V) color. Thus, 
the main sequence in a V, (V—I) CMD is not as tight and 
shows greater intrinsic scatter than the V, (B — V) CMD 
because the unresolved binaries are separated in (V — I) 
color from the single-star main sequence more than they 
are in (B — V) color. An additional selection of the main 
sequence in the V, (V — I) CMD provides a sample of 
probable cluster members with reduced contamination 
by unresolved binaries. For determining the additional 
main sequence selection in the V, (V — I) CMD, we trace 
by eye the blue boundary of the main sequence. The red 
main sequence boundary is defined by offsetting the blue 
boundary by 0.08 mag in color. A 0.15 mag color offset is 
needed to select most of the stars that meet the original 
main sequence selection based on the V, (B — V) CMD 
alone. This additional main sequence selection reduces 
the stellar sample to 584 from the original 1004. The 
best fit isochrone using stars that meet the main sequence 
selection criteria in both CMDs has parameters [Fe/H] = 
-0.05, Age = 1.052 Gyr, (m - M)„ = 12.33 mag, and 
Ay = 0.62 mag. The above solution is shown as the 
empty square in Figure [SJ 

To define the 1-a systematic error on the fitted parame- 
ters, we use the difference between the best fit parameters 
and the parameters determined during the above discus- 
sion of the three sources of systematic error: uncertainty 



in Ry, photometric calibration and binary star contam- 
ination. For the 1-a systematic error in the theoretical 
isochrones, we use the relative difference between our 
isochrone fit parameters for the Hyades and the quoted 
values from the literature. To derive an overall system- 
atic error in the cluster parameters, these four sources 
of systematic error are added in quadrature. The re- 
sulting overall 1-a systematic errors are crpe/H] = 0.08, 
fiAgo = 0.09 Gyr, <7( m _ M ) = 0.12 mag, and a Av = 0.09 
mag. The photometric calibration and Ry uncertainty 
dominate the systematic uncertainty in metallicity. The 
systematic uncertainty in the isochrone dominates the 
systematic uncertainties in the age and distance. Unre- 
solved binary contamination and photometric calibration 
dominate the uncertainty in the reddening. 

3.2. Comparison with Other Determinations 

The best fit metallicity we derive from isochrone fitting 
to NGC 1245 agrees with our independent metallicity 
determination using spectroscopy of individual red-giant 
members (J. Marshall et al., in preparation). Medium- 
resolution spectroscopic indices calibrated from high- 
resolution spectroscopy indicate [Fe/H] = — 0.06± 0.12, 
where the error is the 1-a syste matic error in the sp ectral- 
index-metallicity calibration. iWee fc Lee I (^996) mea- 
sure the metallicity of NGC 1245 using Washington pho- 
tometry and obtain [Fe/H] = —0.04 ± stat 0.05 ± S yst 
0.16. We note that NGC 1245 is commonly quoted 
in the literature as having a super-so l ar me tallicity of 
[Fe/H] = +0.14 as given by the lLvngal <)1987[) open clus- 
ter database 4 . The origin of this high metallicity for 
NGC 1245 is unclear; the source for the high metallicity 
as given in iLvngal l|1987f ) does not contain observations 
for or even discuss NGC 1245. 

Our best fit age for NGC 1245 i s in agreement wit h 
the Age = 1.1 ± 0.1 Gy r fo und bvlWee fc Lee I (11996ft . 
ICarraro fc Patat 1 1)1994 and lSubramaniam I (2003) find 
a youn ger age of Age = 800 My r for NGC 1245. In the 
case of ICarraro fc Patat I l)1994l) . the younger age most 
likely results from their assu mption of [Fe/H] = + 0.14 
for the cluster. In the case of lSuhramaniaml l)2003[) . the 
difference in color between the main sequence turnoff and 
the red-giant clump req uires the younger age. We believe 
the red-giant clump of iSubramaniam (2003) is too red 
by 0.35 mag in (B — V) color (as discussed in section 
Our smaller color difference between the main sequence 
turnoff and the red-giant clump results in the older age 
for NGC 1245. 

The main differences between this study and previ- 
ous investigations of NGC 1245 are that we find a lower 
overall reddening and we find no evidence for significant 
differential reddening. We find E(B -V) = 0.21 ± 0.03, 
where the error includes the systematic uncertainty in 
Ry and Ay. Previous studies find highe r reddenings 
E(B - V) = 0.28 ± 0.03 (IWee fc Lee II1996D . marginally 
consistent with our determination. The previous claims 
for differential reddening across NGC 1 245 have not been 
highly significant and even conflicting. ICarraro fc Patatl 
( 1994) determine a north/south (B — V) color gradient of 
0.04 mag arcmin"" 1 wit h redder colors toward the south, 
and l)Wee fc Lee 1119 96) find a north/south (B — V) color 
gradient of 0.03±0.16 mag arcmin -1 , except the color 

4 http://vizier.u-strasbg.fr 
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gradient is in the opposite sense, with redder colors to 
the north. 

We follo w the same procedu re as in ICarraro fc Patat I 
l|1994[) and I Wee fe Lee I (|1996|) to determine the differen- 
tial reddening across NGC 1245 using the 1.3m dataset. 
To quantify the differential reddening across NGC 1245, 
we start with the same sample of stars that were used 
for the isochrone fit. We further restrict the sample to 
stars with 15.0 < V < 17.0 mag. The differential red- 
dening is modeled as a linear trend between the (B — V) 
color of stars as a function of north/south position. Wc 
formally find a north/south (B — V) color gradient of 
0.0042± 0.0006 mag arcmin -1 with redder colors to the 
north. The color gradient corresponds to a total change 
of A(B — V) ~ 0.03 over 9 arcmin. The color gradient 
in the east/west direction is half the north/south color 
gradient with a similar error. Therefore, we find a color 
gradient that is smaller by an order of magnitude than 
previous investigations. Furthermore, although formally 
significant, we are not convinced of the color gradient's 
reality because the two-dimensional aperture correction 
map and PSF variations both contain errors of order 0.02 
mag. 

4. RADIAL PROFILE 

Stellar encounters drive an open cluster toward 
equipartition of kinetic energ y resulting in mass segrega - 
tion of the clu ster members dBinnev fc Tremainel ll987'). 
The survey of iNilakshi et al. I l)2002j) demonstrates that 
the prominent cores of open clusters generally contain 
the most massive stars, but 75% of a cluster's mass lies 
in a surrounding corona containing only low mass mem- 
bers of the cluster. In this section, we search for the 
signatures of mass segregation in NGC 1245 by study- 
ing its radial surface-density profile as a function of the 
magnitude of the cluster stars. 

The first step in deriving the radial profile is to locate 
the cluster center. Using the stellar positions for stars 
with high-quality photometry shown in the 1.3m CMD 
(Figure [5J, we spatially smooth the number counts with 
a Gaussian of radius ~ 75". We do not use the 2.4m 
dataset to calculate the cluster center because gaps in the 
CCD array bias the derived center. The open circle in 
Figure n denotes the cluster center. Our derived cluster 
center is loca ted 48" east and 75" s outh of the center 
determined bv lCarraro fc Patat 1 1)1994|) a nd 36" east and 
36" so uth of the center determined by iSubramaniam I 
(2003). Since neither study gives exact details on how 
the cluster center was determined, we cannot determine 
the significance of these differences in the cluster center. 
However, we believe that our estimate for the cluster 
center is likely to be more robust, due to our larger field 
of view and fainter limiting magnitude. 

We construct a radial surface-density profile by deter- 
mining the surface density of stars in concentric annuli 
of width ~ 25". The top left panel in Figure H3 shows 
the surface-density profile derived from both the 1.3m 

1/2 

and 2.4m datasets. We approximate the errors as Nf , 
where Ni is the number of stars in each bin. Area in- 
completeness sets in due to the finite size of the detector 
for annuli > 7.3' in the 1.3m dataset and > 11.6' for 
the 2.4m dataset, thereby increasing the errors. Gen- 
erally, the surface-density profiles derived from the two 
datasets agree well, except near the center of the cluster 



and a small overall scaling. The scaling difference results 
from the fainter limiting magnitude in the 2.4m dataset 
(see Figure|5J|. The discrepancy in the inner annuli arises 
from the fact that the cluster center falls very close to 
a gap in the CCD array for the 2.4m detector. There- 
fore, incompleteness affects the first few radial bins of the 
2.4m surface-density profile. We disregard the first two 
radial bins for all surface-density profiles derived from 
the 2.4m dataset. 

As expected, the radial profile of the cluster exhibits 
a gradual decline to a constant field-star surface density. 
We fit the surface-density profile to the model, 



£(r) = E/ + £ 



, v 
1+ I - 

r c 



(4) 



where £/ is the surface density of field stars, So is the 
cluster surface density at r — 0, and r c is the core radius. 
In fitting Equation^] we require a positive S f . We fit this 
profile to both the 1.3m and 2.4m datasets. Although 
the 2.4m surface-density profile extends to larger radii, 
the incompleteness near the center generally results in 
less well-constrained parameters, in particular the core 
radius r c . 

Both fits are excellent: for the 1.3m dataset, x 2 = 15.9 
for 28 degrees of freedom (dof), whereas for the 2.4m 
dataset, x 2 = 39.3 for 36 dof. Table |3] tabulates the 
best fit parameter values and 1-a uncertainties for all 
datasets. A parabola fit to the envelope of x 2 values as 
a function of the parameter of interest determines the 
parameter's error. The solid line in the top left panel of 
FigureElshows the best fit model to the 1.3m profile. The 
logarithmic slope of the surface density is (3 — 2.36±0.71, 
and the core radius is 3.10 ± 0.52 arcmin. The field and 
central surface densities are E/ = 4.0 ±0.8 arcmin -2 and 
So = 14.7 ± 2.9 arcmin -2 , respectively. At the derived 
distance to the cluster, d — 2850 pc, the conversion from 
angular to physical distances is 1.207 arcmin/pc. Thus, 
NGC 1245 has a core radius of ~ 2.6 pc. The core ra- 
dius agrees with the previous deter mination of 2.7± 0.13 
arcmin from [Nilakshi e t al. I (j2002fl . 

To investigate mass segregation, we next divide the 
high-quality photometry sample of stars into two sub- 
samples: bright stars with V < 17.8 and faint stars with 
V > 17.8. There are roughly an equal number of cluster 
stars in each subsamplc. The surface-density profiles for 
both the 1.3m and 2.4m data are shown in the bottom 
left and middle left panels of Figure H3 We find a signifi- 
cant difference between the best fit logarithmic slopes of 
the bright and faint subsamples. For the 1.3m dataset, 
we find (3 = 3.9 ± 1.3 for the bright sample, and we find 
f3 = 1.3±0.9 for the faint sample. This provides evidence, 
at the ~ 2 — a level, that the faint cluster stars in NGC 
1245 are considerably more spatially extended than the 
brighter stars. The fainter stars being spatially extended 
is almost certainly due to the effects of mass segregation, 
which we expect to be important for a cluster of this 
mass and age. We provide additional evidence for mass 
segregation in Section [S] and discuss the implications of 
this mass segregation more thoroughly in Section |SJ 

The radial profile of NGC 1245 implies that the faint 
cluster members may actually extend significantly be- 
yond our field of view. The radial profile for the 1.3m 
data predicts that, at 1-a, 40 to 100% of the faint 
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(V > 17.8) stars outside 10' of the cluster center and 
within the field of view arc cluster members. The strong 
cluster contribution in the outer periphery of our ob- 
served field hampers our ability to obtain an accurate 
statistical subtraction of the field star contamination nec- 
essary to determine the mass function for NGC 1245. We 
discuss this additional complication and our solution in 
Section [5] 

5. MASS FUNCTION 

Typically an assumed mass-luminosity relation allows 
the transformation of the observed lu minosity functio n to 
derive the present day mass function l|Chabrier 12003ft . In 
contrast, we assign a mass to each star individually using 
our best fit isochrone. Minimizing a star's individual xf 
(Equation [2J using the best fit isochrone parameters de- 
termines the star's mass. To determine the mass for each 
star, we use the high quality 2.4m photometry shown in 
the bottom panel CMDs of Figure [5] as it covers a much 
larger field of view. We reduce field star contamination 
of the cluster mass function by restricting the sample to 
stars with xf < 0.04. This xf cutoff selects a region 
in the V, (B — V) CMD roughly similar to the main se- 
quence selection used in the isochrone fit but naturally 
selects stars based on their proximity to the isochrone in 
all three passbands. 

Despite the large field of view of the 2.4m data, the 
radial profile of the cluster suggests that there exists a 
significant number of low mass cluster members out to 
the limits of our field of view. Lacking observations of 
a field offset fr om the cluster, the B esancon theoretical 
galaxy model ijRobin fc Creze I H986) provides an alter- 
native method to correct the observed mass function for 
field star contamination. We obtain from the Besangon 
galaxy model electronic database the stellar properties 
and photometry for a one-square-degree field centered on 
the cluster's galactic coordinates without observational 
errors. A spline fit to the median photometric error in 
0.5 mag bins in each passband allows us to simulate the 
observed errors in the theoretical galaxy field photome- 
try. We add a zero-mean, unit-standard-deviation Gaus- 
sian random deviate scaled to the observed photometric 
error model to the theoretical galaxy field photometry. 
We apply the same error selection criteria (a < 0.08) as 
was applied to the observed CMD and place a saturation 
cutoff at V=14.0 mag to the theoretical CMD. 

To make a qualitative comparison between the theo- 
retical field CMD and the observed CMD, we adjust the 
theoretical CMD for the smaller observed 2.4m field of 
view. In the comparison, a star is included in the the- 
oretical CMD only if a uniform random deviate is less 
than the ratio between the two field of views. Figured 
compares the CMDs of the theoretical galaxy model in 
the top panels and the observations for stars beyond 12.7 
arcmin of the cluster center for the 2.4m data in the bot- 
tom panels. The dashed line shows the best fit isochrone 
for reference. As shown in Figure the theoretical and 
observed field CMDs are qualitatively very similar. 

To calculate the field star contamination in the mass 
function, we first obtain the mass for the stars in the 
theoretical field CMD and apply the same xf < 0.04 se- 
lection in the same manner as the observed field CMD. 
For both the theoretical field CMD and observed cluster 
CMD, we calculate mass functions using even log-spaced 



intervals of 0.071 log(M©). Number counts from the en- 
tire one-square-degree field of view provide the basis for 
the theoretical galaxy mass function. The larger field of 
view of the theoretical galaxy model than the observed 
field requires scaling of the number counts and corre- 
sponding Poisson error of the theoretical galaxy mass 
function by the ratio of field of views. Subtracting the 
scaled theoretical galaxy field number counts from the 
observed number counts provides the cluster mass func- 
tion. The cluster mass function is shown as the heavy 
solid line in the upper panel of Figure |S1 The top axis 
shows the median apparent V magnitude for the cor- 
responding mass bin. The vertical dashed line delin- 
eates mass bins that are complete to better than 90%. 
Comparing the cluster mass function to a mass function 
calculated by subtracting the theoretical galaxy number 
counts without applying the photometric error selection 
criteria determines the completeness of a mass bin. The 
long-dashed line labeled 'Salpeter' in Figure illustrates 
the typical observed mass function slope, a = —1.35. 
A power law in linear-mass, linear-number-count space 
fit with weights to the mass function bins not affected 
by incompleteness yields a slope steeper than Salpeter, 
a = -3.12 ±0.27. 

To verify the validity of using the theoretical galaxy 
field in place of an observed off-cluster field, we cal- 
culate an alternative mass function using the observed 
field outside 12.7 arcmin of the cluster center for field 
contamination subtraction. The entire 2.4m field of 
view again provides the basis for the observed number 
counts. The smaller field of view of the outer periphery 
requires scaling the field contamination number counts 
and corresponding Poisson error by the relative area fac- 
tor, / = 6.4. The light-solid line in the upper panel of 
Figure |H1 shows the resulting cluster mass function using 
the outer periphery for determining the field contami- 
nation. The mass bin positions are offset slightly for 
clarity. The two mass functions agree within 1 — a for 
all bins with log M > but disagree significantly at the 
low mass end. The theoretical galaxy subtracted mass 
function contains greater number counts at the low mass 
end, consistent with our conclusion that the outer re- 
gions of our observed field of view contain a significant 
number of low mass cluster stars. Both mass functions 
contain a drop in the number counts at log(M) = 0.071. 
We question the significance of this drop in the number 
counts. The drop occurs at a magnitude with some of 
the heaviest field star contamination, making the field 
star contamination subtraction less certain. 

To further the evidence for mass segregation in NGC 
1245, we calculate the mass function inside and outside 
a radius of 4.2 arcmin from the cluster center. According 
to the cluster radial profile, this radius roughly separates 
half of the cluster members. The inner and outer mass 
functions are shown as the solid and short-dashed lines, 
respectively, in the lower panel of Figure^ The fit to the 
inner region mass function yields a shallower slope that is 
consistent with Salpeter, a = — 1.37± 0.2. Ignoring the 
highest mass bin, the best fit slope for the inner region 
is even shallower, a = —0.56 ± 0.28. We find that the 
outer region is highly enriched in low mass stars. The 
outer region mass function has a very steep slope, a = 
— 7.1± 1.2. The mass bins with the down pointing arrows 
have negative values, and the point represents the 2-cr 
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upper confidence limit. Since the mass function is fit in 
linear space, the fit includes bins with negative values. 

The radial profile and mass function demonstrate that 
a significant fraction of the low mass cluster members 
reside in the outer periphery of the cluster. Previous 
studies of NGC 1245 generally find Salpeter and shal 
lower slopes for the mas s function iCarraro fe Patat 



ll994t|Subramaniam l2003f) . As warned in lSubramamam 
( 2003), the smaller field of view of these previous studies 
from which to determine the field contamination biases 
the mass function against low mass members that reside 
in the outer periphery of the cluster. The mass func- 
tions from the previous studies agree well with our mass 
function for the interior of NGC 1245, but we conclude 
that a steeper slope than Salpeter, a — —3.12 ± 0.27, is 
more appropriate for NGC 1245 as a whole down to the 
completeness limit of M = 0.8M Q . 

Our deeper CMD obtained from combining the MDM 
2.4m data (FigureOJ) allows us to extend the mass func- 
tion to even lower masses. Determining the mass func- 
tion to lower masses improves our mass estimate of the 
cluster, calculated in the following section. To calculate 
the mass function, we follow the same procedure as out- 
lined above. At these fainter magnitudes, the theoretical 
galaxy field overpredicts the number of stars. Thus we 
must resort to using the outer periphery of the MDM 
2.4m field of view for estimating the field contamination. 
The heavy solid line in Figure shows the resulting mass 
function based on the deeper CMD. The light solid line 
reproduces the equivalent mass function from Figure |H1 
based on the shallower CMD data. The vertical dashed 
line represents the completeness limit at M — 0.56M Q . 
The deep mass function has a slope, a = —0.5 ± 0.4, sig- 
nificantly steeper than the slope for the shallow mass 
function, a = —2.8 ± 0.8. The significantly differing 
slopes suggest there is a turnover in the mass function 
at log(M/M Q ) = -0.05. Fitting to the last three com- 
plete mass bins in the deep mass function yields a slope, 
a = 1.8 ± 0.5. The low mass end of the deep mass func- 
tion likely does not fall as steeply as measured, since we 
use the outer periphery of the MDM 2.4m field of view 
to quantify the field contamination, which likely contains 
many low mass cluster members. Thus, the measured 
slope at the low mass end represents an upper limit. 

6. TOTAL MASS AND DYNAMICS 

Numerous dynamical processes take place in open clus- 
ters. Interactions between the cluster members dramat- 
ically affect binary and planetary systems, mass segre- 
gation with in the cluster, and the timescale for cluster 
disso l ution ijHu rley fc Ijhara |l2Q Q2albl : iGiersz fc Heggiel 
119971 iPortegies Zwart et al. 1120011) . The cluster's ini- 
tial total mass sets the timescale for these dynamical 
processes and governs its dynamical evolution. Unfortu- 
nately, several difficulties prevent a completely empirical 
determination of the total mass of NGC 1245. First, 
as we show in Section the surface-density profile sug- 
gests the cluster extends significantly beyond our field of 
view. Second, a substantial population of low mass clus- 
ter members exists below the completeness limit of our 
mass function. Third, we do not resolve binary systems. 
Therefore, the total mass in observed cluster members 
(~ 820M Q ) only places a lower limit to the current clus- 
ter mass. Finally, even overcoming the heretofore diffi- 



culties, stellar evolution and evaporation results in ad- 
ditional complications for determining the initial birth 
mass of the cluster. 

In order to provide a rough estimate of the current total 
cluster mass, we must adopt a number of assumptions in 
order to extrapolate the observed mass to the total mass. 
We must account for stars outside our survey area and 
stars fainter than our magnitude limit. The total cluster 
mass determined by integrating the radial profile model 
we fit in Section 0] slowly converges and is unrealistically 
large for an integration to infinity. In reality, the tidal 
field of the Galaxy truncates the outer radius of the clus- 
ter. Therefor e, we r efit the observed dataset to the King 
profile l)King 111962ft . 
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1+1^ 



-1/2' 



(5) 

where rt is the tidal radius of the cluster. Note that, as 
r t — * oo, this profile reduces to the radial profile given in 
Equation^lfbr (3 = 2. For a given set of parameters; Eo, 
r c , and rt; the integral over Equation determines the 
total number of stars, 



v 2v j'ln(l+z t )-[3(l + zt) 
iv = 7rr c 2_,o 



1/2 _ 



1][(1 



xt 



,1/2 _ 



l+X t 

(6) 

where x t = (r t /r c ) 2 l|King l H.962'). The total mass is then 
M = (m} b s N, where {m) b s = LOOM© is the average 
mass of observed cluster stars. 

Unfortunately, King profile fits to either dataset yield 
only lower limits for the value of r t . The lower limit 
only for r t is primarily due to the covariance between rt 
and £y. Therefore, we need an additional constraint on 
either the surface density of field stars or the tidal radius 
in order to determine the total cluster mass. Attempting 
to constrain the fit using values of E/ determined from 
the theoretical galaxy field star counts (see Section 
results in unrealistically large values of r t > 40 arcmin. 
Thus, we adopt a different approach. 

The tidal radius of the cluster is set by both the local 
tidal field of the Galaxy and the total mass of the cluster. 
We can estimat e the tidal radius as the radius of the 
Lagrange point l)Binnev fc Tremainelll987f ). 



n 



(—) 

\3MmwJ 



1/3 



D, 



(7) 



where D is the distance of the cluster from the center of 
the Galaxy, M is the total mass of the cluster, and Mmw 
is the mass of the Galaxy interior to D. For our derived 
distance of d — 2850 pc, assuming R = 8 kpc, and 
Galactic coordinates (I — 147°.6,6 = — 8°.9), the Galac- 
tocentric distance of NGC 1245 is D = 10.50 kpc. For a 
flat rotation curve with v — 220 km s , we find Mmw — 
1.18 x 10 n M Q , and thus r t ~ 20 pc(Af/2500M Q ) 1 / 3 . 

We use equations © and 10 to solve for the total 
cluster mass M and tidal radius r t while simultaneously 
constraining r c ,£^, and Eo from fitting Equation [S] to 
the observed radial surface-density profile. However, we 
must first account for low mass stars that are below our 
detection limit. To do this, we use the best fit mass 
function (dN/dlogM oc M a , with a = -3.12) to the 
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turnover in the mass function at O.85M . For stars be- 
tween the hydrogen-burning limit and 0.85M©, we as- 
sume a = 1.0. Thus, our assumed cluster mass function 
is 

^ / A/" 3 ' 12 , M > O.85M , 
dlogM \ Af 10 , O.O8M < M < O.85M K > 

The adopted mass function slope at the low mass end, 
a = 1.0 is a compromise between our measured lower 
limit to the slope, a = 1.8 (see Section |2J, and other 
clusters that have an observed mass function slope as 
shallow as a ~ 0.4 at the low mass end jP risinz ano et al. I 
I200H iBouvier et aTll2003j) . This yields a total cluster 
mass (in our field of view) of ~ 1312M and a correction 
factor of 1.61 to the observed mass in cluster stars. 

We use an iterative procedure to determine M and rt- 
We assume a value of r t and then fit the King profile 
(Equation 03) to the observed radial density profile. The 
best fit King profile yields trial values of r c , £/, Soj and 
the total number of cluster stars via Equation|S| We then 
apply the correction factor to account for stars below our 
completeness limit to determine the total cluster mass, 
M — 1.61(m) ob s N '. We then use Equation [7| to predict 
a new value for r t . The procedure iterates until conver- 
gence with a tidal radius of r t ~ 20 arcmin (16.5 pc) 
and a total mass of Af = 1312 ± 9OAf , where the error 
is the 1-cr statistical error from the fit. We determined 
the total cluster mass using the shallow 2.4m dataset as 
this is the dataset used to determine the mass correction 
factor and surface-density profile. 

The mass correction factor, sensitive to the mass func- 
tion slope at the low mass end, dominates the uncer- 
tainty in the cluster mass. For example, varying the log- 
arithmic slope of the mass function below O.85M in the 
range 0.6 < a < 1.4 changes the derived total mass by 
+2OO-14OAf . 

Our derived mass does not account for the presence of 
unresolved binaries. Without knowledge of the proper- 
ties of the binary population, it is difficult to assess the 
magnitude of the correction. However, for a large binary 
fraction, the correction could be as large as 50%. We 
have also not attempted to correct the derived present- 
day mass for evolutionary effects, such as evaporation 
or stellar evolution, to derive an initial cluster mass. 
For our adopted mass function (Equation |SJ, the cor- 
rection to the total cluster mass due to stellar evolu- 
tion is small, ~ 12%. About 10% of the clus ter stars 
are lost to the Galaxy per relax ation time l)Spitzer I 
119871 iPortegies Zwart et al. 1 120011) . Since the current 
age of NGC 1245 is approximately 8 times larger than 
the current relaxation time (see below), it is possible 
that the cluster has lost an appreciable amount of mass 
over its lifetime. The numer ical N-body simulations of 
IPortegies Zwart et al. I i|2001|) predict the initial cluster 
mass based on the current number of red-giant stars. Us- 
ing their Equation 1, an age of 1 Gyr for NGC 1245 and a 
total of 40 observed red-giant stars yields an initial total 
cluster mass, M D = 3120Af Q . 

The relevant quantity that determines the timescale for 
dynamical evolution of a stellar clu ster is the half-mass 
relaxation timescale (jSpitzer Ill987|) . 

/ r 3 \ 1/2 N l/2 



Here Th m is the radius enclosing half of the total mass of 
the cluster, (m) ~ O.63M is the average mass of cluster 
members, and In A ~ In 0.4iV is the Coulomb logarithm. 
This is roughly the average time over which the velocity 
of a typical cluster member changes by order unity due 
to random encounters with other cluster members. For 
NGC 1245, r hm = 3.8 pc, and t rlx ~ 130 Myr. Since the 
age of NGC 1245 is ~ 1 Gyr, we expect this cluster to 
be dynamically quite evolved. Therefore, the evidence of 
mass segregation found in Sections 0] and [S] is not sur- 
prising. 

7. CONCLUSION 

In this paper, we report our results based on BVI 
photometry of the open cluster NGC 1245 using the 
MDM 1.3m and 2.4m telescopes. Here, we improve 
on the observations for this cluster by covering a six- 
times greater area, obtaining the first CCD /-band pho- 
tometric data, and acquiring the first CCD data un- 
der photometric conditions. With these significant im- 
provements over earlier studies, we are able to deter- 
mine more precisely the physical parameters of the clus- 
ter as a whole, rather than just the inner cluster core. 
Our photometric data resolve some of the confusion in 
the literature regardi ng the correct photometric zero- 
point for this cluster ijWee fc Lee 1119961: iSubramaniam I 
2003). B ased on isochron e fits employing the Y 2 cal- 
culations <Yi et al. ll200l|) . we confirm the findings of 
iWee fc Lee I l)1996|) and (J. Marshall et al., m prepara- 
tion) that this cluster has a slightly subsolar metallicity, 
[Fe/H] = -0.05 ± 0.03 (statistical) ± 0.08 (systematic). 
Our best fit age is 1.04 ± 0.02 ± 0.09 Gyr. In contrast to 
previous studies, we do not find evidence for significant 
differential reddening. We find a U-band extinction of 
A v = 0.68 ± 0.02 ± 0.09, and with the aid of 2MASS K s 
photometry, we constrain the ratio of absolute to selec- 
tive extinction to be Ry = 3.2 ±0.2. The resulting abso- 
lute distance modulus is (m- Af) = 12.27±0.02±0.12. 

With the large field of view provided by the MDM 
8K Mosaic imager, we find NGC 1245 is a highly re- 
laxed cluster; a majority of the low mass cluster mem- 
bers reside in an extended halo. The mass function slope, 
a = -3.12 ± 0.27, down to M = O.85M is steeper than 
the Salpeter value of a — —1.35, found by earlier stud- 
ies. The previous studies did not have a sufficient field of 
view to detect the low mass cluster members that pref- 
erentially reside in the outer periphery of the cluster. 
Based on the observed stellar surface-density profile and 
an extrapolated mass function, we derive a total cluster 
mass, M = 1300 ± 90 ± 17OM . 

This paper is the first in a series for the Search 
for Transiting Extrasolar Planets in Stellar Systems 
(STEPSS) project. The main goal of the STEPSS 
project is to determine the fraction of stars with close- 
in, Jupiter-sized planets using the photometric transit 
technique. By searching for extrasolar planets in open 
clusters of known metallicity, stellar density, and age, we 
hope to constrain the impact of these parameters on the 
resulting planetary systems. The accurate physical pa- 
rameters for NGC 1245 derived in this paper will be used 
in the analysis of the 19 night run targeting NGC 1245 
to search for transiting extrasolar planets (B. S. Gaudi et 
al., in preparation). The NGC 1245 transit search, con- 
sisting of 5 minute time resolution J-band photometry 
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Date (2001) #Exps Filter Exp (s) FWHM (arcsec) Comments 



Oct. 27 


9 


V 


300 


1.4 


overcast 




9 


I 


300 


1.1 




Nov. 7 


7 


V 


300 


1.8 


partly cloudy 




7 


I 


300 


1.2 




8 


5 


B 


240 


1.6 


cirrus 




5 


V 


240 


1.5 






5 


I 


240 


1.3 




10 


5 


B 


240 


1.4 


cirrus 




5 


V 


240 


1.6 






5 


I 


240 


1.3 




11 


5 


B 


120 


1.1 


clear 




5 


V 


120 


1.1 






5 


I 


120 


1.0 






1 


1 


240 


0.9 





taken over 19 nights, will also provide valuable informa- 
tion on the variable star content of NGC 1245 (J. Pepper 
et al., in preparation). 

This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and Space 



Administration and the National Science Foundation. 
This work was supported in part by NASA through a 
Hubble Fellowship grant from the Space Telescope Sci- 
ence Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under 
NASA contract NAS5-26555; the Menzel Fellowship from 
the Harvard College Observatory; and by NASA grant 
NAG5-13129. 
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Fig. 1. — A false color image of the open cluster NGC 1245 obtained with the MDM 8K Mosaic imager (the 4x2 array of CCDs are 
numbered) on the MDM 2.4m telescope. The large inscribed square is the field of view for the MDM 1.3m telescope data obtained with 
a single CCD imager. The blue, green, and red color channels of the false color image consist of combined images in the BVI passbands, 
respectively. The circle shows the cluster center. 



Table 2. MDM 1.3m Observations 



Date (2002) #Exps Filter Exp (s) FWHM (arcsec) Comments 



Jan. 31 


2 


B 


120 


2.0 


non-photometric 




3 


V 


60 


1.9 






3 


1 


30 


1.8 






3 


1 


20 


1.8 




Feb. 5 


6 


B 


120 


1.7 


non-photometric 




4 


V 


60 


1.6 






6 


1 


20 


2.0 




Feb. 7 


2 


B 


90 


1.6 


photometric 




2 


V 


60 


1.1 






3 


1 


20 


1.7 
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Fig. 2. — Color magnitude diagrams (CMDs) of the open cluster NGC 1245. Error bars represent the median sample standard deviation 
in the stellar color as a function of magnitude. Top panels: Data from the MDM 1.3m. Bottom panels: Data from the MDM 2.4m. Left 
panels: V versus (B — V). Right panels: V versus (V — I). 



Table 3. Fitted Parameters and I — a Errors to Surface-Density Pro- 
files. 



Datasct 


Sample 




So 


Tc 





T,f (model) 






# / arcmin 


# / arcmin 2 


arcmin 


# / arcmin 2 








(#/pc 2 ) 


(#/pc 2 ) 


(pc) 


(#/pc 2 ) 


1.3m 


All 


4.05 ± 0.84 


14.72 ± 2.88 


3.10 ±0.52 


2.36 ± 0.71 


5.087 ±0.040 






(5.90 ± 1.22) 


(21.43 ±4.20) 


(2.57 ±0.43) 




(7.401 ± 0.058) 


1.3m 


Bright 


2.29 ±0.23 


7.61 ± 1.52 


3.27 ±0.51 


3.87 ± 1.34 


2.189 ±0.024 






(3.33 ± 0.34) 


(11.08 ± 2.21) 


(2.72 ± 0.42) 




(3.187 ±0.035) 


1.3m 


Faint 


0.49 ± 1.93 


7.76 ±3.47 


4.61 ± 2.42 


1.29 ± 0.88 


2.887 ±0.028 






(0.71 ± 2.80) 


(11.29 ±5.05) 


(3.82 ± 2.00) 




(4.203 ± 0.041) 


2.4m 


All 


4.70 ± 1.10 


18.89 ± 12.19 


2.22 ±1.44 


1.66 ± 0.78 


5.440 ± 0.039 






(6.84 ± 1.61) 


(27.51 ± 17.8) 


(1.84 ± 1.19) 




(7.920 ± 0.057) 


2.4m 


Bright 


1.59 ±0.28 


7.58 ±3.03 


2.75 ±0.95 


2.53 ± 1.06 


2.053 ±0.025 






(2.32 ± 0.40) 


(11.04 ±4.41) 


(2.28 ± 0.79) 




(2.989 ± 0.036) 


2.4m 


Faint 


1.54 ±3.14 


29.17 ±43.21 


0.40 ±5.01 


0.70 ± 1.29 


3.379 ± 0.030 






(2.24 ± 4.57) 


(42.46 ± 62.91) 


(0.33 ±4.15) 




(4.919 ± 0.044) 
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Fig. 4. — The heavy solid line represents the best fit isochrone solution overlying the color magnitude diagrams based on the MDM 1.3m 
data. The light solid lines show the main sequence selection boundaries. 
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Fig. 5. — Contour lines represent the 1-, 2-, & 3-sigma confidence regions on the joint variation in metallicity and age. The filled square 
shows the best fit solution. The open square shows the best fit solution after reducing binary star contamination. The open triangle shows 
the best fit solution assuming a ratio of selective-to-total extinction of Ry = 3.0 instead of Ry = 3.2. The starred point shows the best fit 
solution assuming an /-band calibration fainter by 0.02 mag. 
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Fig. 6. — Left panels show the radial surface-density profiles of stars in the NGC 1245 field for three different samples of stars. The top 
row shows the radial profile for all stars. The middle panel shows the profile for the subset of these stars with V < 17.8, while the lower 
panel is for V > 17.8. The solid circles are derived from the 1.3m data, while the open squares are from the 2.4m data. In each panel, the 
solid line shows the best fit model to the 1.3m data, which is composed of a cluster profile (dotted line) plus a constant field-star surface 
density (dashed line). The panels on the right show the 68%, 95%, and 99% confidence regions on the fitted logarithmic slope of the density 
profile j3 and core radius for the 1.3m data. The filled point shows the best fit model for the 1.3m data, whereas the open square is the 
best fit for the 2.4m dataset. The shaded region is the 68% confidence region for the fit to the 2.4m dataset. In the upper-right panel, we 
also show the 68% confidence regions for the bright and faint samples. 
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Fig. 7.— Top Pan els: Besangon theoretical galaxy model V, (B - V) and V, (V - I) CMDs for the Galaxy position of NGC 1245 
IRobin fc C rczc 198t|. Photometric errors and selection criteria are included in the theoretical galaxy to match the observed CMDs shown 
in the bottom panels. Bottom Panels: Observed V, (B — V) and V, (V — I) CMDs for stars beyond 12.7 arcmin of the cluster center using 
the 2.4m data. In all panels, the dashed line shows the best fit isochrone. 
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Fig. 8. — Top Panel: Mass function using a theoretical galaxy model (heavy solid line) and outer periphery of the field of view (light solid 
line) for determining the field contamination. The vertical line delineates mass bins complete to better than 90%. Bottom Panel: Mass 
function for stars inside (solid line) and outside (dashed line) 4.3 arcmin. The down pointing arrows represent mass bins with negative 
values, and the point on the tail represents the 2-sigma upper limit. 
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Fig. 9. — Heavy solid line shows the mass function calculated from the deep 2.4m data (Figure^ using the outer periphery of the field 
of view for determining the field contamination. The vertical line delineates mass bins complete to better than 90%. The thin solid line 
reproduces the equivalent mass function based on the shallow 2.4m data as shown as the dashed line in the upper panel of Figure Isl 



